Pareto Law in a Kinetic Model of Market with Random Saving Propensity 
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We have numerically simulated the ideal-gas models of trading markets, where each agent is 
identified with a gas molecule and each trading as an elastic or money-conserving two-body colli- 
sion. Unlike in the ideal gas, we introduce (quenched) saving propensity of the agents, distributed 
widely between the agents (0 < A < 1). The system remarkably self-organizes to a critical Pareto 
distribution of money P(m) ~ m~^ +1 ^ with v ~ 1. We analyse the robustness (universality) of 
the distribution in the model. We also argue that although the fractional saving ingredient is a 
bit unnatural one in the context of gas models, our model is the simplest so far, showing self- 
organized criticality, and combines two century-old distributions: Gibbs (1901) and Pareto (1897) 
distributions. 
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Considerable investigations have already been made to 
study the nature of income or wealth distributions in 
various economic communities, in particular, in differ- 
ent countries. For more than a hundred years, it is 
known that the probability distribution P(m) for income 
or wealth of the individuals in the market decreases with 
the wealth m following a power law, known as Pareto law 
Q: 

P(m) oc m-( 1+u \ (1) 

where the value of the exponent v is found to lie be- 
tween 1 and 2 0, 0, Q • It is also known that typically 
less than 10% of the population in any country possesses 
about 40% of the wealth and follow the above power law. 
The rest of the low-income group population, in fact the 
majority, clearly follows a different law, identified very 
recently to be the Gibbs distribution IE • Studies 
on real data show that the high income group indeed fol- 
low Pareto law, with v varying from 1.6 for USA [||, to 
1.8 — 2.2 in Japan [3j. The value of v thus seem to vary 
a little from economy to economy. 

We have studied here numerically a gas model of a 
trading market. We have considered the effect of saving 
propensity of the traders. The saving propensity is as- 
sumed to have a randomness. Our observations indicate 
that Gibbs and Pareto distributions fall in the same cat- 
egory and can appear naturally in the century-old and 
well-established kinetic theory of gas : Gibbs distribu- 
tion for no saving and Pareto distribution for agents with 
quenched random saving propensity. Our model study 
also indicates the appearance of self-organized criticality 
in the simplest model so far, namely in the kinetic 
theory of gas models, when the stability effect of savings 
[To| is incorporated. 

We consider an ideal-gas model of a closed economic 
system where total money M and total number of agents 



N is fixed. No production or migration occurs and the 
only economic activity is confined to trading. Each agent 
i, individual or corporate, possess money m,(i) at time 
t. In any trading, a pair of traders i and j randomly 
exchange their money [E 0, [n|, such that their total 
money is (locally) conserved and none end up with neg- 
ative money (rrii(t) > 0, i.e, debt not allowed): 

mi(t) + mj[t) = m»(i + 1) + mj(t + 1); (2) 

time (t) changes by one unit after each trading. The 
steady-state (t — ► oo) distribution of money is Gibbs one: 

P(m) = (l/T)exp(-m/T);T = M/N. (3) 

Hence, no matter how uniform or justified the initial 
distribution is, the eventual steady state corresponds to 
Gibbs distribution where most of the people have got 
very little money. This follows from the conservation of 
money and additivity of entropy: 

P(m 1 )P(m 2 ) = P{mi + m 2 ). (4) 

This steady state result is quite robust and realistic too! 
In fact, several variations of the trading, and of the 'lat- 
tice' (on which the agents can be put and each agent 
trade with its 'lattice neighbors' only), whether compact, 
fractal or small- world like 0] , leaves the distribution un- 
changed. Some other variations like random sharing of 
an amount 2m2 only (not of mi + 7712) when mi > m 2 
(trading at the level of lower economic class in the trade), 
lead to even drastic situation: all the money in the mar- 
ket drifts to one agent and the rest become truely pauper 

nana- 

In any trading, savings come naturally 10]. A saving 
propensity factor A is therefore introduced in the same 
model (see for model without savings), where each 
trader at time t saves a fraction A of its money m; (t) and 
trades randomly with the rest: 
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rrii(t+l) = mi(t) + Am; m,j{t + l) = m,j(t) — Am (5) 
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where 

Am = (1 - X)[e{Tm(t) + m 3 {t)} - mi {t% (6) 

e being a random fraction, coming from the stochastic 
nature of the trading. 

The market (non- interacting at A = and 1) becomes 
'interacting' for any non- vanishing A(< 1): For fixed A 
(same for all agents) , the steady state distribution Pf (m) 
of money is exponentially decaying on both sides with 
the most-probable money per agent shifting away from 
m = (for A = 0) to M/N as A -> 1 (Fig. 1(a)). 
This self-organizing feature of the market, induced by 
sheer self-interest of saving by each agent without any 
global perspective, is quite significant as the fraction of 
paupers decrease with saving fraction A and most peo- 
ple end up with some fraction of the average money in 
the market (for A — » 1, the socialists' dream is achieved 
with just people's self-interest of saving!). Interestingly, 
self-organisation also occurs in such market models when 
there is restriction in the commodity market |14|. Al- 
though this fixed saving propensity does not give yet the 
Pareto-like power-law distribution, the Markovian na- 
ture of the scattering or trading processes (eqn. |@J|) 
is lost and the system becomes co-operative. Indirectly 
through A, the agents get to know (start interacting with) 
each other and the system co-operatively self-organises 
towards a most-probable distribution (m p ^ 0). 




FIG.l: Steady state money distribution (a) Pf(m) for the fixed 
(same for all agents) A model, and (b) Pj(rn) for some typical 
values of A in the distributed A model. The data is collected from 
the ensembles with N = 200 agents. The inset in (b) shows the 
scaling behavior of Pf(m). For all gents start with average 

money per agent M/N = 1. 

In a real society or economy, A is a very inhomogeneous 
parameter: the interest of saving varies from person to 
person. We move a step closer to the real situation where 
saving factor A is widely distributed within the popula- 
tion. One again follows the same trading rules as before, 
except that 

Am = e(l - \j)mj{t) - (1 - Ai)(l - e)m l (t) (7) 

here; A^ and Xj being the saving propensities of agents i 
and j. The agents have fixed (over time) saving propensi- 
ties, distributed independently, randomly and uniformly 



(white) within an interval to 1 (see [15] for preliminary 
results): agent i saves a random fraction Aj (0 < A, < 1) 
and this Xi value is quenched for each agent (Xi are in- 
dependent of trading or t). Starting with an arbitrary 
initial (uniform or random) distribution of money among 
the agents, the market evolves with the tradings. At each 
time, two agents are randomly selected and the money 
exchange among them occurs, following the above men- 
tioned scheme. We check for the steady state, by looking 
at the stability of the money distribution in successive 
Monte Carlo steps t. Eventually, after a typical relax- 
ation time (~ 10 5 for N = 200 and uniformly distributed 
A) dependent on N and the distribution of A, the money 
distribution becomes stationary. After this, we average 
the money distribution over ~ 10 3 time steps. Finally 
we take configurational average over ~ 10 5 realizations 
of the A distribution to get the money distribution P(m). 
It is found to follow a strict power-law decay. This de- 
cay fits to Pareto law (|TJ) with v = 1.02 ± 0.02 (Fig. 2). 
Note, for finite size N of the market, the distribution has 
a narrow initial growth upto a most-probable value m p 
after which it falls off with a power-law tail for several 
decades. As can be seen from the inset of Fig. 2, this 
Pareto law (with v ~ 1) covers the entire range in m of 
the distribution P(m) in the limit N — » oo. We checked 
that this power law is extremely robust: apart from the 
uniform A distribution used in the simulations in Fig. 2, 
we also checked the results for a distribution 

p(X) ~ |A - A| a , A ^l, 0<A<1, (8) 

of quenched A values among the agents. The Pareto law 
with v = 1 is universal for all a. The data in Fig. 2 
corresponds to Ao = 0, a — 0. For negative a values, 
however, we get an initial (small m) Gibbs-like decay in 
P(m) (see Fig. 3). 
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FIG. 2: Steady state money distribution P(m) in the model for 
distributed A (0 < A < 1) for N = 1000 agents. Inset shows that 
the most probable peak m p shifts towards (indicating the same 
power law for the entire range of m) as N — * oo: results for four 
typical system sizes N = 100, 200, 500, 1000 are shown. For all 
cases, agents play with average money per agent M/N = 1. 




FIG. 3: Steady state money distribution P(m) in the model for 
N = 100 agents with A distributed as p(X) ~ A" with different 
values of a. The inset shows the region of validity of Pareto law 
with 1 + v = 2. For all cases, agents play with average money per 
agent M/N = 1. 
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FIG. 4: Money of the i-th trader with time. For fixed A for all 
agents in the market: (a) with A = 0, (b) A = 0.5, (c) A = 0.9. 
For distributed A (0 < A < 1): for agents with (d) A = 0.1, (c) 
A = 0.5 and (f) A = 0.9 in the market. All data are for N = 200 
(M/N = 1). 



In case of uniform distribution of saving propensity A 
(0 < A < 1), the individual money distribution Pf(m) 
for agents with any particular A value, although differs 
considerably, remains non-monotonic: similar to that for 
fixed A market with m p (X) shifting with A (see Fig. 1). 
Few subtle points may be noted though: while for fixed 
A the rrip(A) were all less than of the order of unity (Fig. 
1(a)), for distributed A case m p (A) can be considerably 
larger and can approach to the order of N for large A 
(see Fig. 1(b)). The other important difference is in the 
scaling behavior of P/(m), as shown in the inset of Fig. 
1(b). In the distributed A ensemble, Pf(m) appears to 
have a very simple scaling: 



P f (m) ~ (1 - A).F(m(l - A)), (9) 




for A — > 1, where the scaling function F(x) has non- 
monotonic variation in x. The fixed (same for all agents) 
A income distribution Pf(m) do not have any such com- 
parative scaling property. It may be noted that a small 
difference exists between the ensembles considered in Fig 
1(a) and 1(b): while f mPf(m)dm = M (independent 
of A), J mPf(m)dm is not a constant and infact ap- 
proaches to order of M as A — > 1. There is also a marked 
qualitative difference in fluctuations (see Fig. 4): while 
for fixed A, the fluctuations in time (around the most- 
probable value) in the individuals' money m,i(t) gradually 
decreases with increasing A, for quenched distribution of 
A, the trend gets reversed (see Fig. 4). 



FIG. 5: Money distribution in cases where saving propensity A is 
distributed uniformly within a range of values: (a) When width 
of A distribution is 0.5, money distribution shows power-law for 
0.5 < A < 1.0; (b) When width of A distribution is 0.2, money 
distribution starts showing power-law when 0.7 < A < 0.9. Note, 
the exponent u ~ 1 in all these cases when the power law is seen. 
All data are for N = 100 (M/N = 1). 

We now investigate on the range of distribution of the 
saving propensities in a certain interval a < Aj < 6, 
where, < a < b < 1. For uniform distribution within 
the range, we observe the appearance of the same power 
law in the distribution but for a narrower region. As may 
be seen from Fig. 5, as a — > b, the power-law behavior is 
seen for values a or & approaching more and more towards 
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unity: For the same width of the interval \b — a\, one gets 
power-law (with same v) when 6 — > 1. This indicates, 
for fixed A, A = corresponds to Gibbs distribution, and 
one gets Pareto law when A has got non-zero width of 
its distribution extending upto A = 1. This of course 
indicates a crucial role of these high saving propensity 
agents: the power law behavior is truely valid upto the 
asymptotic limit if A = 1 is included. Indeed, had we 
assumed An = 1 in ffl. the Pareto exponent v immedi- 
ately switches over to v = 1 + a. Of course, A / 1 
in (JSJ leads to the universality of the Pareto distribu- 
tion with v = 1 (independent of An and a). Indeed this 
can be easily rationalised from the scaling behavior JJjJ: 
P(m) ~ Jo Pf(m)p(\)d\ ~ m- 2 for p(X) given by © 
and m~( 2+ ") if An = 1 in (JSJ (for large m values). 
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FIG. 6: Cumulative distribution Q(m) = P(m)dm of wealth 
m in USA 0] in 1997 and Japan 3] in 2000. Low-income group 
follow Gibbs law (shaded region) and the rest (about 5%) of the 
rich population follow Pareto law. The inset shows the cumulative 
distribution for a model market where the saving propensity of the 
agents is distributed following JSJ with An = and a = —0.7. The 
dotted line (for large m values) corresponds to v = 1.0. 

These model income distributions P(m) compare very 
well with the wealth distributions of various countries: 
Data suggests Gibbs like distribution in the low-income 
range Q (more than 90% of thepopulation) and Pareto- 
like in the high-income range |3| (less than 10% of the 
population) of various countries (Fig. 6). In fact, we 
have compared one model simulation of the market with 
saving propensity of the agents distributed following JSJl , 
with Ao = and a = —0.7. This model result is shown 
in the inset of Fig. 6. The qualitative resemblance of the 
model income distribution with the real data for Japan 



and USA in recent years is quite intriguing. In fact, for 
negative a values in ©, the density of traders with low 
saving propensity is higher and since A = ensemble 
yields Gibbs-like income distribution @, we see an ini- 
tial Gibbs-like distribution which crosses over to Pareto 
distribution with v = 1.0 for large m values. The po- 
sition of the crossover point depends on the magnitude of 
a. The important point to note is that any distribution 
of A near A = 1, of finite width, eventually gives Pareto 
law for large m limit. The same kind of crossover behav- 
ior (from Gibbs to Pareto) can also be reproduced in a 
model market of mixed agents where A = for a finite 
fraction of population and A is distributed uniformly over 
a finite range near A = 1 for the rest of the population. 

We also considered annealed randomness in the saving 
propensity A: here Aj for any agent i changes from one 
value to another within the range < Aj < 1, after each 
trading. Numerical studies for this annealed model did 
not show any power law behavior for P(m); rather it 
again becomes exponentially decaying on both sides of a 
most-probable value. 

We have numerically simulated here ideal-gas like mod- 
els of trading markets, where each agent is identified with 
a gas molecule and each trading as an elastic or money- 
conserving two-body collision. Unlike in the ideal gas, 
we introduce (quenched) saving propensity of the agents, 
distributed widely between the agents (0 < A < 1). For 
quenched random variation of A among the agents the 
system remarkably self-organizes to a critical Pareto dis- 
tribution (JTJ of money with v ~ 1.0 (Fig. 2). The ex- 
ponent is quite robust: for savings distribution p(X) ~ 
An — A |", A 7^ 1, one gets the same Pareto law with 
v = 1 (independent of Ao or a). It may be noted that 
the trading market model we have talked about here has 
got some apparent limitations. The stochastic nature of 
trading assumed here in the trading market, through the 
random fraction e in © , is of course not very straightfor- 
ward as agents apparently go for trading with some def- 
inite purpose (utility maximization of both money and 
commodity). We are however, looking only at the money 
transactions between the traders. In this sense, the in- 
come distribution we study here essentially corresponds 
to 'paper money', and not the 'real wealth'. However, 
even taking money and commodity together, one can ar- 
gue (see [lij) for the same stochastic nature of the trad- 
ings, due to the absence of 'just pricing' and the effects 
of bargains in the market. 

Apart from the intriguing observation that Gibbs 
(1901) and Pareto (1897) distributions fall in the same 
category and can appear naturally in the century-old and 
well-established kinetic theory of gas, that this model 
study indicates the appearance of self-organized critical- 
ity in the simplest (gas) model so far, when the stability 
effect of savings incorporated, is remarkable. 
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